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Abstract 


We derive a closed expression for the spatial shift experienced by a black soliton colliding with a 
shallow dark soliton in the context of deformed non-linear Schrodinger models. A perturbative scheme is 
developed based on the expansion parameter 1 / (•yv) « 1 , where v is the velocity of the incoming shallow 
dark soliton, 7 = 1/^/1 — u 2 /uf and v s is the Bogoliubov sound speed, therefore it is not restricted to 
small deformations of the integrable NLS model. As applications of our formalism we consider the 
integrable NLS model and the non-integrable cubic-quintic NLS model with non-vanishing boundary 
conditions. Extensive numerical simulations are performed in order to verify our results. A variant of 
the analysis for gray-gray soliton collision is discussed regarding a fast broad soliton and a slow thin 
soliton. 



1 Introduction 


The mathematical theory of soliton collisions in the integrable models is a well developed subject. In partic¬ 
ular, the solitons emerge with their properties completely unchanged after collision between them, except for 
relevant phase shifts. However, certain non-linear field theory models possess solitary wave solutions and it 
is difficult to know a priori if they are in fact true solitons with analytical multisoliton solutions. Regarding 
this issue the so-called quasi-integrability concept has recently been put forward related to the deformations 
of the integrable sine-Gordon (SG) and non-linear Schrodinger (NLS) models [1]. According to this new 
concept certain deformed integrable field theory models share some properties with their truly integrable 
counterparts for some configurations, for example they possess an infinite number of quasi-conserved charges. 
For the case of two solitary wave collision those charges are asymptotically conserved in the scattering pro¬ 
cess, although the charges vary in time their values in the far past and the far future remain the same. The 
both numerical and analytical tools were used to define and describe the concept of quasi-integrability. 

Our aim is to predict the results of solitary wave collisions in deformed NLS models with non-vanishing 
boundary conditions. Such theories appear in diverse areas of non-linear science, such as condensed matter 
physics, plasma physics and, in particular, in the study of Bose-Einstein condensates. We are also motivated 
by the above mentioned quasi-integrability property in deformed NLS models with non-vanishing boundary 
conditions in which sufficiently fast solitary waves effectively pass through each other. 

In the context of analytical calculations some results have been obtained only for the cases of small 
perturbations of the NLS model [2,3,4]. For nearly integrable models it has been implemented a perturbation 
theory allied to the the inverse scattering transform (1ST) method [5]. These methods have been applied 
mainly to the NLS solitons that decay at infinity, i.e. the bright solitons. They encompass the above 
mentioned perturbation theory based on 1ST and a direct perturbation theory [6]. On the other hand, the 
dark solitons are intensity dips on a continuous wave background with a phase change across their intensity 
minimum. However, the non-vanishing boundary of dark solitons, which may change when the deformations 
are present, introduces serious complications when applying the perturbative techniques developed for bright 
solitons. In the last years, various improvements have been made to the calculations and methods [7, 8, 9] 
based on the so-called complete set of squared Jost solution (eigenvectors of the linearized NLS operator) in 
order to apply to the non-vanishing boundary conditions and dark soliton perturbations. The implementation 
of direct methods, however, consider small perturbations around the NLS model [4, 10]. Moreover, in [4] it 
has been argued that the squared Jost functions associated with the dark soliton might be an insufficient 
basis, so rendering problematic the issue of a complete basis for the solution space of the deformed NLS 
equation. 

Here we propose a method to study shallow solitary wave collisions in the context of a perturbation 
theory performed order by order in an expansion parameter related to the both velocity of the incoming 
soliton and the Bogoliubov sound speed. The expansion parameter has a kinematic origin and does not 
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appear explicitly in the Lagrangian of the model. In this way our series expansion is not restricted to small 
deformations around the integrable NLS model. We will assume that at zeroth order in this parameter 
the colliding solitary waves effectively pass through each other and provide the leading order contribution 
to the spatial shift experienced by a stationary soliton after colliding with an incoming soliton. Similar 
perturbative framework has recently been considered in the study of scattering of ultra-relativistic kink-type 
solitary waves of the deformed SG model [11, 12]. We present the relevant developments in the context of a 
non-Lorentz invariant field theory associated to the NLS complex field and show that the absolute value of 
the spatial shift experienced by a black soliton after collision with an shallow gray soliton decreases as the 
velocity increases, in this way resembling the behavior of ultra-relativistic kink-(anti)kink collisions. 

It has recently been shown that the perturbed dark soliton of the NLS model comprises an inner region, 
containing the core of the soliton and a moving shelf, and an outer region which matches to the boundary 
at infinity [13, 4]. The main features of those results relevant to our discussions here are the appearance 
of a shelf of finite length around the both sides of the soliton the edge of which propagates with speed 
determined by the background intensity, and that the difference between the amplitude of the shelf and the 
outer background is very small compared with the amplitude of the background intensity. 

The main characteristics of our perturbative scheme are the following. We consider a linear superposition 
of two solitary wave solutions as a good background solution during collision, such that the interaction of 
the colliding solitons occurs significantly in the overlapping region which is a region contained in the inner 
region of the perturbed stationary solitary wave. The continuous wave background will be subtracted on 
top of the linear superposition in order to maintain the non vanishing boundary condition at infinity. So, 
we are interested in investigating the effects of short-range interactions between the dark solitons and the 
long range contributions arising from the shelves mentioned above will be ignored in the integrals employed 
to compute the spatial shift. In the rest frame of the stationary solitary wave where the other approaches 
with velocity |i>| —> v s , the space-time integration of the overlapping region becomes proportional to 
where 7 = 1/y/l — ( v/v s ) 2 >> 1, v s being the Bogoliubov sound speed. This property is systematically 
used to provide a perturbative framework in order to calculate the spatial shift experienced by the stationary 
soliton, using the small value of (^-) as the parameter of our perturbation expansion. The perturbation 
parameter (^) arises after space-time and suitable field space transformations are performed simultaneously. 
In this regard the implementation of the perturbative analyses of the non-relativistic deformed NLS model 
is slightly different from the relativistic deformed SG model in which an analogous factor appears directly 
from a coordinate transformation to the rest frame of the static soliton [ 11 ]. 

The paper is organized as follows. In section 2 we discuss a deformed NLS model with non-vanishing 
boundary conditions and the relevant equations leading to dark soliton solutions. In section 3 we develop a 
perturbative scheme in order to compute the spatial shift experienced by a black soliton after colliding with 
a gray soliton. In section 4 we present some examples of application of our formalism. The integrable NLS 
model is considered in some detail in order to check our perturbative scheme. Afterwards, the collision of 
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dark solitons of the cubic-quintic NLS model is studied as an example of a non-integrable theory. In section 
5 we develop numerical simulations for the both models considered in section 4. In section 6 we discuss the 
main results and possible future directions of research. Finally, in the appendix A the properties of dark 
soliton interaction in the integrable NLS model are presented, and in the appendix B some properties of 
dark solitons of the cubic-quintic NLS model are discussed 


2 Deformed NLS models and solitary waves 


Various generalizations and deformations of the integrable NLS model arise in physical applications. We will 
consider the deformed NLS models of the type 


d Id 2 


df[\^(x,t)\ 2 } 

d\ip(x,t)\ 2 


ip(x,t) = 0, 


( 2 . 1 ) 


where / : 1R + —> ]R. The model (2.1) allows dark soliton type solutions in analytical form for certain 
functions /[/]; for example, the integrable NLS case corresponds to /[/] = /3/ 2 /2 (see Appendix A) and 
the non-integrable cubic-quintic NLS (CQNLS) defined by /[/] = j3I 2 /2 — a/ 3 /6 [14] (see Appendix B). 
Among the models with saturable non-linearities [15], the case /[/] = ^ p s (I + ) also exhibits analytical 

dark solitons [16]. The qualitative properties of traveling waves of the NLS model for general non-linearities 
and non-vanishing boundary conditions have been studied in [17]. The formalism developed here to study 
soliton collisions will require the deformed NLS model to possess dark-soliton type solutions obtained either 
analytically or numerically, as well as the condition that the linear superposition of two such solitons provides 
a good background solution during their collision. 

We will consider translationally invariant solutions of the form 


ijj{x, t) = cj)(z)exp[id(z) + iwt], z = x — vt, (2.2) 

satisfying the following non-vanishing boundary condition ( nvbc) 


\${,Z t ihoc)| — IV’ol, dz$(z— >-±oo) — 0, d z 9( z — ^± 00 ) — 0* 


(2.3) 


The nvbc (2.3) is suitable for wave solutions in the form of localized ‘dark’ pulses created on the continuous 
wave background ( cwb ) at rest 

Mx,t) = \He i(wt+S) , (2.4) 

where w = — / , [|V’o| 2 ] (here and in the upcoming sections f'[y\ = is the phase shift experienced by 

the cwb with intensity IV’ol 2 , in the non-linear medium with intensity \i/j\ 2 and varying according to the law 
~ d ^d[t\ 2 ^ • Moreover, we will consider a family of systems such that small oscillations around |^>| = |i/’o| of 
(2.1) defines a hyperbolic system with sound speed 

Vs = \MVf"[\M 2 ], n\H 2 ] > 0 . ( 2 . 5 ) 
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2.1 Solitary waves 

Among the solutions of type (2.2) let us assume some dark solitons as follow 


M*,t) = M^)e iwt e iC , (2.6) 

V> 2 (s,i) = Mz)e ie(z) e iwt , z = x - vt, (2.7) 


where %l>i will be the stationary ‘black’ soliton and if) 2 will be the incoming ‘gray’ soliton. <f> 1, 4> 2 ,9 are 
smooth functions of their arguments such that </>i : 1R —>• [^IV’ol, |V'oI] 1 02 ■ H —> [?0)1"001] C 1R + , 9 : 1R —>• 
[9~ ,9 + ], Moreover, they converge exponentially fast to the non-vanishing boundary conditions: 0i(±oo) —> 
A | if>o |, 02(±oo) —> |-0015 0(±oo) —> 9 ± . The constant phase ( in (2.6) has been introduced for later purpose, 
and the frequency w has been assumed to be the same as the one of the cwb in (2.4). 

One can write an equation for <j> 1 by substituting the Ansatz for 0i (2.6) into (2.1) 

0' 1 '-2( W + /'[0?])0 1 =O. (2.8) 


From the equation (2.8) one can get the relationships 


d(j>i 

dx 


±V2 

±F[cf> 1 


'IV’ol 2 


(fV}-f'Uo\ 2 ])dI 


.1 V2 


(2.9) 


The corresponding solution for the NLS model is a black soliton 0i, the solution (A.4) with v = 0, 
which is an antisymmetric function of space and develops an abrupt n phase shift and zero intensity at its 
center. The =F signs in (2.9) correspond to the regions x € [—00, Xo] and x £ [a;o,+oo], respectively, where 
Xq represents the zero intensity position of the black soliton. The relevant black soliton <f> 1 of the CQNLS 
model is presented in (B.ll) which presents quite similar properties. 

Substituting the Ansatz for ip 2 (2.7) into (2.1) one gets the following equations 

02-2 ( w + f(4\)fo + v 2 (fa~ 1 -^) = 0 ( 2 . 10 ) 

±{ti[9'(z)-v}} = 0 , ( 2 . 11 ) 

where w = —/'[IV’ol 2 ]- Integrating the last equations one gets 


d(j> 2 
dz 


d9_ 

dz 


9(x, t) 


±V2\ {f[i]~fm 2 ])di- 0 

Ibol 2 1 

±G[0 2] 

H[<t> 2] 

, r 02 m 

Lo\ G[y] y 
±0[02] 


v 2 (0l - |0o| 2 ) 2 


1/2 


( 2 . 12 ) 


(2.13) 


(2.14) 
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In the NLS case one has the solution (A.4) such that <j >2 = \i(>(x,t)\ and 9 being the space dependent 
phase of ip(x,t). Whereas, the gray soliton (f> 1 of the CQNLS model is presented in (B.4) and the space 
dependent phase in (B.6). Notice that in the both models the =F signs in (2.12) correspond respectively to 
the regions z € [— 00, z 0 ] and z € [ zq , +00] of the dark soliton, where zo is the coordinate in which the soliton 
amplitude approaches its minimum intensity (dip). 

The analytical expression for the two-dark soliton solution of the integrable NLS model is presented in 
(A.5) which allows the computation of the exact value of the spatial shift developed after collision of black 
and gray solitons (A. 10). In Fig. 1 we have plotted a collision process of these solitons. 

The above relationships (2.9) and (2.12)-(2.14) will be used below in our study of the soliton collisions in 
the framework of a perturbative scheme and in the Appendix B in order to find the CQNLS dark solitons. 

3 Perturbative expansion for black and gray soliton collision 

In order to study dark soliton collisions we will assume that the system (2.1) possesses some solitary waves 
such that their linear combination provides a good background solution during collision. For variety of 
functionals /[/] and in the relevant ones used in the examples below this property is verified. Consider a 
solution before, during and after the collision of two-dark solitons 

^(x,t) = ipi(x,t) + ipz(x,t) + h(x,t ) - ipo(x,t). (3.1) 

Here ijji(x,t) represents the stationary ‘black’ soliton (2.6) and ipz(x,t) the moving ‘gray’ soliton (2.7). 
The function h(x, t) describes the perturbation effects around the core of the stationary black soliton gen¬ 
erated by the collision. It has been included the continuous wave background (2.4) in order to subtract the 
nvbc reference term ipo(x, t) on top of the superposition, since the functional /[/] associated to the potential 
of the system is not periodic. The constant phase parameters ( and 5 associated to the black soliton (2.6) 
and the cwb (2.4), respectively, will play below an important role in order to compute the spatial shift. The 
Fig. 2 shows a gray-black soliton system of the CQNLS model before, during and after collision. Notice 
the successive deformations experienced by the static soliton during the entire process, in particular the Fig. 
2.c. shows the deformation, described by the function h(x,t) above, of the core of the static soliton during 
the nearest position of the both solitons. 

As mentioned above a dark soliton under perturbation develops a shelf the edge of which propagates out 
on the both sides of the soliton. In order to study the radiation shed by an evolving dark soliton one can 
linearize the NLS equation (2.1) about the (cwb) as ip a (x,t) « \4>o\^ wt + \ipOa\e xwt ■ So, the moving shelves 
contributions associated to each soliton can be incorporated as ip a (x, t) + \ipo a \e lwt (a = 1,2); however it has 
been shown that IV’Oal << If/’ol [13, 4], therefore we will neglect them in our computations of the spatial shift 
at first order in l/(i>7). 
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The equation (2.1) can be linearized in order to find the equation for h(x,t). So, one has 

id t h + — Wi(x,t)h — W%(x, t)h — W$(x,t)h = S(x,t) 

+ AWi(x, t)h + Albjfi, t)h + AW|(x, t)h (3.2) 

where 


Wi(x,t) 

= rm\ 

(3.3) 

W^{x,t) 

= f"m 2 WA 2 , 

(3.4) 

W%(x,t) 


(3.5) 

S(x, t) 

= (/'[l*o| 2 ] - /'[|^i| 2 ])^i + (/'Pol 2 ] - f'[ifa])'h 



P/'Po| 2 P/'Po| 2 ]Po 

(3.6) 

HO 

< 

= nM-rm 2 }, 

(3.7) 

AWi(x,t) 

= /"[|^o| 2 ]|^o| 2 -/"[|^i| 2 ]|V>i| 2 , 

(3.8) 

A W£{x,t) 

= /"Po| 2 ]Po) 2 -/"[|Vh| 2 ]V^ 

(3.9) 

^0 (x,t) 

= i>l{x,t) + 1 p 2 (x,t) - 'lpo(x,t). 

(3.10) 


The expressions W\ , IT2 1 and W'j in the l.li.s. of (3.2) are the potential terms for the perturbation around an 
isolated stationary dark soliton -f/q; and AWi, ATI ^ 1 and AW.f are their relevant changes due to the incoming 
dark soliton ip2- We refer to S(x, t) as the external source, whose real and imaginary parts behave as step-like 
functions around the regions occupied by the solitons, becoming a single step-like function when the two 
solitons overlap as it was shown in the Fig. 3 for the collision of gray and black solitons of the CQNLS 
model. Notice that in the absence of the gray soliton one can make h = 0, ip2 —>■ V’o, then S, AWi, AW^ 
and AW.2 vanish identically, whereas they are non vanishing in the regions occupied by the gray soliton if) 2. 
Considering the incoming dark soliton as the source of perturbation one notices that h, S, AWi , A W \, AW2 
are non vanishing when the two dark solitons overlap. 

There is another equation for h, where the bar stands for complex conjugation, which can be obtained 
by taking the complex conjugate of (3.2). 

In order to find solutions of the eq. (3.2) we expand h(x, t ) as 

h(x, t) = 9x(t)u\(x) (3-11) 

A 

where {ma(*)} form an orthonormal basis and solve the eigenvalue equation in matrix notation 



(3.12) 


-R 2 + / , [|^i| 2 ] + / ,, [|^| 2 ]|^i 

/"[ IVfil 2 ]«h) 2 


/"[ l^flWh) 2 

-3^ + / , [IV’i| 2 ] + /"[|V’i| 2 ]|V’i 
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Figure 1: (color online) The initial condition (5.5)-(5.6) for a dark-black soliton 
system of the NLS model (4.1) is plotted for /3 = 4, |-0o I = 3, £ = 1.15026 as a 
continuous line (green). The gray soliton initially located at X 2 = — 8 travels to 
the right with velocity v = 4.64. The black soliton is located at X\ = 4.2. The 
final configuration after collision is the dotted line (red). Note the spatial shift 
Ax ~ —0.124 experienced by the black soliton. 

A remarkable fact is that one can find the first zero modes (A = 0) 

ul(x) = e l %ipi, ul(x) = d x ipi, (3-14) 

where ipi solves the stationary version of the deformed NLS equation (2.1) for any /[/]. These zero modes are 
associated to the phase and translational symmetries of the model, respectively. Regarding this issue, let us 
mention that the linearized eigenvalue equation around a dark soliton solution in the integrable NLS model 
possesses four linearly independent zero modes u J 0 (x,t),j = 1,4 and their exact analytical expressions 
have been presented in [18]. Those zero-modes are associated to the following four fundamental symmetries 
of the NLS model: phase, translational, Galilean and scale invariance, respectively. 

3.1 spatial shift and a perturbative expansion 

Next, we make use of the linearized equation (3.2) in order compute the spatial shift experienced by the 
black soliton after collision with a gray soliton. It is considered that the shallow gray soliton is approaching 
the static black soliton from the negative x direction as in Figs.l and 2. Let us expand the stationary black 
soliton as 


ip i{x + A x(t),t) = ipi(x,t) + d x 'ipi(x,t)Ax(t) + ... 
Then in order to extract the spatial shift we take 

h(x,t ) = d x ipi(x,t)Ax(t) + ... 


(3.15) 


(3.16) 
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Figure 2: (color online) a) Before collision. The initial red profiles of the CQNLS dark solitons. The left 
soliton is the gray soliton and the right one is the static black soliton. b) The gray soliton is approaching 
the black soliton from the left, c) The nearest position of the both solitons. For comparison the initial 
profile of the black soliton is plotted as dotted green line, d) The gray soliton pulse is starting to emerge 
to the right from the position of nearest collision. The initial black soliton is the green curve, e) After 
collision. The final profile of the dark solitons. The right soliton is the gray soliton and the left one is 
the static black soliton. f) The spatial shift of the black soliton. The final soliton (red curve left) has a 
position shifted with respect to the initial black soliton position (green curve right). The parameters are 
a = 0.6, p = 12.0, IV'ol = 2.0, v = 6.07, v s = 6.19677. 
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Figure 3: (color online) The real (continuous curves) and imaginary (dashed curves) parts of the 
source S(x,t) generated by the incoming gray soliton through its interaction with the black soliton 
of the CQNLS model. The initial profile (5.7)-(5.8) has been used for the gray-black soliton system. 
The gray soliton approaches the black soliton from the negative direction of the coordinate x. a) 
The initial profile of the source. The left step-like functions are around the gray soliton position 
and the right ones are around the static black soliton. b) The source S during the collision and the 
nearest position of the both solitons. c) The source profile after the collision. 
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The function d x ipi(x,t) above is precisely the zero mode Uq(x) (3.14) associated to the space translation 
symmetry of the deformed NLS model. Comparing (3.16) with the eigenmode expansion as in (3.11) one 
can see that the ellipses represent the expansion terms in the eigenmodes which are orthogonal to d x ipi(x). 
Regarding the velocity change of the black soliton AC one can consider a term like h(x, t) = Uq(x)AV + ..., 
where Uq(x) is the zero mode associated to Galilean symmetry of the deformed NLS model. However, as 
we will show below the spatial shift does not present a time dependence at leading order of perturbation 
theory, so the velocity change AC must be zero at this order. This fact is also verified through our numerical 
simulations for some concrete examples. We are mainly interested in the effects of the zero mode associated 
to space translation and postpone for a future work the treatment of the remaining zero inodes and higher 
order effects. 

Since AlCi, AW] and AW| are multiplied by h in (3.2) their effect are small compared to the source 
term S, so substituting (3.16) into (3.2) and considering only the source term S one can get 

id x ipi (x, t) = S(x,t) (3.17) 

Projecting out the last relation and its complex conjugate onto the zero modes associated to space 
translation symmetry and <9 x '0i, respectively, and integrating with respect to x one gets 


dAx(t) 

dt 


1 


m 

/•+oo 


r+oo 


dx 


d x ipiS(x,t) - d x ij)\S{x,t) 


V(t) = 2 / dx d x ipid x i/>i. 

J — OO 

Consider the initial condition Ax{t — oo) = 0, then integrating (3.18) w.r.t. t one has 

A7(f) 


A x(t) = 


m 


Af(t) = i dr 


r +oo 


dx 


’ —OO «/ —OO 


d x ipiS(x,T) - d x ipi S(x,t) 


(3.18) 

(3.19) 


(3.20) 

(3.21) 


The last equation, taking into account the source S (3.6) and the expressions of the dark solitons (2.6)-(2.7), 
can be rewritten as 


m 

l*o| 2 


/ t r+oo 

dr / dx 
-oo J —OO 


#if 

dx l L 


<i> 2 sin (9 — C) + IV'ol sin (C + <5)J /'[|^o| 2 ] - 
(f >2 sin (9 - Qf[\<h\ 2 ] - IV’ol sin (( + 6)f'[\ip 0 \ 2 }}, 


+ 24 > 1 (f >2 cos (6 - C) - 20i |^)o | cos (C + 5)- 20 2 IV^o I cos (9 - 8). 


(3.22) 

(3.23) 


In the next steps we will change the space-time (x,t) integration to the field space (</>i,</> 2 ) integration 
in the limit t —> +oo. So, consider the coordinate transformation 

x' = x, z = (x — vt), dxdr = —dx'dz (3.24) 


10 



The above transformation and (2.9)-(2.12) allow us to write dz = Adfa/Glfa], where the minus (plus) sign 
corresponds to the left (right) interval of integration dcj> 2... (/^^ d<f> 2 •••)> respectively. Notice that qo 

is the minimum intensity of the gray soliton, and it is assumed that the gray soliton is approaching the 
localized region of the black soliton from the negative x direction with a speed v. Therefore, in the limit 
t —> +00 the expressions for Af(+oo) and V(+oo) in field space can be written as 


Af 


V 


4 
v . 


/•Ibol rlbol ^ 

-Ibol U G l^\ 


| (j )2 sin (0 — C) + \ % l ) o\ sin (C + $) f'[\^o\ 2 ] - 


fo sin (6 — C)/ , [|</ , 2 | 2 ] - |V’o|sin(C + <5)/ , [|'!/’o| 2 ]}, 

r-lbol 




-Ibol 


dcf>iF[(f>i], 


(3.25) 

(3.26) 


where the field 6 appearing above explicitly, as well as through the expression for |^o| given in (3.23), must 
be substituted in terms of the ‘integration variable’ <f >2 through 6 = Q[(j) 2 ] from (2.14). 

So, after the colliding shallow ‘gray’ soliton ip 2 is sufficiently far from the stationary ‘black’ soliton t/q 
the spatial shift of it becomes 


Ax = 


Af 
2 M' 



(3.27) 


where Ad = \ f dx[(^-) 2 + (F[4 >^]) 2 ] is the energy of the stationary black soliton. 

For some models which allow analytic solitary wave solutions one observes that the maximum |-i/»o | and 
minimum qo intensities of the incoming wave are associated, respectively, to some roots of the functional 
G[<j> 2 ]. This feature is present in the models with saturable nonlinearities [16] and in the NLS and CQNLS 
models presented in the appendices A and B, respectively. Then, the integrand in (3.25) diverges in the both 
end points of the integration interval, requiring a certain regularization procedure in order to extract a finite 
value for the spatial shift expression. Moreover, the difference between the roots behave as (|'0o| — qo) 0 
for v —> v 3 , i.e. in the sound speed limit. Making use of the above properties (discussed in the appendices 
for the NLS and CQNLS cases) and in order to see the appearance of the perturbative parameter l/(uy) in 
the expression for Af in (3.25) one makes a transformation to a convenient field space 


dd>id(f >2 = 7 — r&[y]d(bidu, 

{vv 


with (f >2 = $[y], y = uj 7 , 7 1 = ^j\ — (y-) 2 , such that (3.25) can be written as 


Af = - 


(try) 

4 

(try) 


/ I’M 

d(f>iZ[<j>i,u], 

-Ibol 
duK, [u] 


where 


(3.28) 


(3.29) 

(3.30) 


Z[ct> 


i) 1 




$[-] sin (0 - C) + IV’ol sin (C + 8) f[ |^ 0 | 2 ] - 
7 J 
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(3.31) 


sin (9 - C)/'[(<^]) 2 ] - \M sin (C + 8)f'[ |^ 0 | 2 ]} 

|*o| 2 = + (4>[-]) 2 + |V>o| 2 + 20!$[-] cos (0-0- 2^|^o| cos (C + 5)- 

7 7 

7/ 

2$[—]|^o| cos (0 — 5). (3.32) 

7 

The functional /C[it] is obtained from (3.29) by integrating the functional Z[<j)i,u] in (j)\. The phase 
function 6 which appears in the integration (3.30) through (3.31)-(3.32) must be written in terms of u by 
making use of the relationship 9 = ©[$[-]] defined in (2.14), where the minus sign has already been taken 
into account in the relevant region of integration. Notice that the factor 1/ (vj) above is not directly obtained 
from a space time transformation as in the analog computations in relativistic scalar field theory models 
[11, 12], in our case it also requires an additional field space transformation (3.28). The parameter u a in 
general may depend on v and the parameters of the theory. For example in the NLS case below we will 
find that u 0 is a vanishing parameter in the sound speed limit and in the CQNLS case the relevant u 0 is a 
constant equal to unity. So, making use of eq. (3.27) and computing the zero’th order term of the integral 
in (3.30) one can write the spatial shift as 

2 r f u ° i 

= d “' C l"l] (0) + e|(»7) -2 ], (3-33) 

where A4 is the energy of the black soliton and [...](o) means that only the zero’th order term of the full 
integral is considered. Remarkably, in order to compute Ax we need only the form of the functional /[/] 
which enters into the relationships (3.29)-(3.33) explicitly and through the functional G\<j> 2 ] as defined in 
(2.12). Despite the fact that the collision is dissipative the spatial shift does not exhibit a time dependence 
implying that no velocity change is observed at leading order. 

In this way the above result for the spatial shift furnishes the leading order term of the stationary soliton 
perturbation h(x,t) through the expression (3.16). Note that in that expression this contribution accounts 
only for the zero mode related to space translation. 

In the construction above the continuous wave background is assumed to be at rest, i.e. k = 0 in (2.4) 

and it is suitable for treating the shallow dark (gray)-black soliton collisions. A variant of the above method 

applied to high speed gray-gray soliton collisions requires a non-rest cwb , as we will see below. We consider a 

small amplitude (shallow) dark soliton with velocity V 2 and a wide width A- colliding with a large amplitude 

gray soliton with velocity (—tq) and a narrow width A-, such that i >2 > v\ (i>i,i >2 > 0) and P 2 < Pi- Notice 

that in the integrable NLS case the analytic 2-dark soliton solution involves these parameters satisfying 

2 _ 2 

the relationship 2k = V 2 — v\ + , where k is the wave number of the background (A.2). In order to 

compute the spatial shift experienced by the slow thin soliton after colliding with the fast broad soliton in 
the context of our formalism it is convenient to consider the co-moving coordinate of the low velocity thin 
soliton. In this frame, the large amplitude soliton is stationary, whereas the small amplitude soliton travels 
with a relative velocity, and the oscillatory cwb is flowing with certain phase velocity. In this frame consider 
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the cwb solution of the deformed NLS model (2.1) as a non-rest oscillatory phase 

Mx,t) = \Me i{kx+wt+S \ (3-34) 

where k, w = — \k 2 — / , [|V , o| 2 ] and w/k are the wave number, the dispersion relation and the phase velocity 
of the background, respectively. A stationary soliton can be transformed into a nonstationary form by using 
the Galilean symmetry. In fact, assuming that ijj(x,t) is a solution of the deformed NLS equation (2.1) and 
c a constant, then tp'(x,t) = e l ^ cx ^ £ ^ t ^(x — ct,t ) will also be a solution. So, consider a solution before, 
during and after the collision of gray-gray solitons in the co-moving frame of the slow soliton as 

W(x,t) = ipi(x,t) + ip 2 (x,t) + h(x,t) - ipo(x,t) 

= ( j) 1 (x)e Wl{x) e i{kx+wt+Cl) + (j) 2 {x - e i(kx+wt+Cs) + h{x,t) 

iv» 0 | y kx+wt+s ). 

Here ipi(x,t) represents the stationary gray soliton and ip 2 (x,t) the moving gray soliton such that v is the 
relative velocity between the soliton and the background. The constant phase parameters £i, ( 2 ,5 have 
been included. The function h(x,t) describes the perturbation effects around the core of the stationary 
gray soliton generated by the collision in the co-moving frame. As above it has been included the cwb 
(3.34) in order to subtract the nvbc reference term ipo(x,t) on top of the superposition. Following similar 
steps as in section 2.1 one can obtain the relationships between the functions 4>\{x) and 9\(x), as well as 
between (f> 2 (z) and 0 2 (z), (z = x — vt), valid in the co-moving frame. The expansion in (3.35) can be 
substituted into the deformed NLS equation (2.1) in order to get a linearized equation for h(x,t) as in 
eqs. (3.2)-(3.10). The perturbative computations must be performed in the co-moving frame and they can 
follow the relevant steps as above; the stationary gray soliton can be Taylor expanded and related to the 
function h(x,t ) as in eqs. (3.15) and (3.16), respectively. Analogous expressions to the ones in eqs. (3.17)- 
(3.23) can directly be reproduced. The relevant space-time integration as in (3.22) would be transformed 
to the field space integration: (x,t) —> {4>i,4> 2 ) in the limit t —> +00. A convenient transformation in 
the field space such that dfyidfo ~ j^d^dcj)^ can be implemented as in (3.28). However, in the above 
construction some integrals would require a renormalization procedure in order to remove the divergences 
associated to the contribution of the oscillatory cwb [15]. For example, the integral (3.19) becomes in this 
case V = 2 dx{[-^(j> i(x )] 2 + }, which is divergent for k =/= 0 since < j ) 1 —> \ipo\ as \x\ +00. 

The renormalized expression V r = 2 dx{[-^(j)i{x)] 2 + fc 2 |V’o| 4 [^yy — j^p]} must be considered. An 

analogous expression to (3.33) will provide the leading order (fq ) -1 contribution to the spatial shift. 

Next, we present some examples of application of the above formalism to the gray-black soliton collisions, 
paying special attention to the relevant field space transformation in each model in order to extract the 
expansion parameter 1 /(^ 7 ). 


(3.35) 

(3.36) 
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4 Examples 


4.1 NLS model 


Let us consider the integrable defocusing NLS model in order to test the formalism developed above. One 
has the NLS model defined by 


- P\i>(x,i)\ 2 i/)(x,t) = 0, 


(4.1) 


where we have considered /'[/] = fdl in (2.1). Note that in this model a closed analytic expression is 
available for the spatial shift of colliding gray-black solitons (A.10)- (A.11). Next let us compute the spatial 
shift applying our formalism. So, take 


m = § i 2 - 


Then, from the eqs. (2.9)-(2.14) one can get 

F[<f>i] = V^d^ol 2 ~ <fi), G[<f>2\ = 100,2 ~ ^ m 2 ] =v(l~y?f- 


9 = ± arcsin ( -4=— ^ — i 


(4.2) 

(4.3) 

(4.4) 


where 9\ is a constant of integration. The function G[p] possesses two roots, i.e. {|p|, v/y/fi}. The 
first root is the maximum intensity, whereas the second one is the minimum intensity (dip) of the dark 
soliton, respectively. In the field space integration of A f (3.25) the field (j>i appears only in the expression of 
/'Pol 2 ] = /3|4> 0 | 2 , where |4/ 0 | 2 is given in (3.23); so after integrating in </>i one has 


M=- 


4|p|/3 


Ibol 


d<t> 2 


w J ib G M 


| ^> 2 sin (6» — C) + |V , o|sin(C + 5) 


4|V’o| 5 


2|pp cos {6- 5) - $2 sin (9 -Q~ Mo| 3 sin (( + <5)j, 


(4.5) 


where the relationship (4.4) must be used in order to write the integrand only in terms of the ‘integration 
variable’ <f> 2 . The above expression can be written in terms of the perturbation parameter (-A) as in (3.29) 
in order to organize the perturbative series. In fact, the expression (4.5) can be written as 


a r = - 


4u s 

(vj) 


du 


'o “o 


l If. — U‘ 


■| <j> 2 [u] sin (0 — C) + |P| sin (£ + 6) 


4|^o| 2 ,2r 1 

- i- <j> 2 M — 


21PPM cos (9 - (5) - 4>l[u] sin (9 - () - |p| 3 sin(C + <5)j, 


PM = 1/-T+7 2 


IMo |/3 


(4.6) 

(4.7) 


Notice that the upper limit of the integration interval u 0 approaches zero in the sound speed limit v —> v s . So, 
due to the overall factor (^L) it can be sufficient to evaluate the zeroth order term of the integral. However, 
it is instructive to perform the full integration in order to see the nature of the divergences. In fact, the 
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integrand in (4.5) presents certain divergences at the both ends of the interval [-4^, |'0o|]- They are associated 
to the roots of the functional G[(f> 2 ] in the denominator, so let us examine the relevant divergent terms by 
performing the full integration. The associated indefinite integral in (4.5), after integration, becomes 


|6wcos (20i + £ + 5) + \/~P I'f/’ol 4 sin (£ + Oi) + 3sin(0i + £ + 2(5) - 3sin(0i - C) } x 
logd^ol 2 - 4>l) + 

UuAlVbl cos(C + Oi) - 3m\//3|V , o| cos (0i - C) + 3vA|^| cos(C + 0i + 25) + 


3 ^\Jv 2 s — v 2 L 

3/3|V'o | 2 sin (C — <5) + 4/3|?/5 0 | 2 sin (C + (5) - 6 w 2 sin (20! + ( + 5) + 3/3\ip 0 \ 2 sin (20! + ( + (5) 


arctanh 




V v s - 


\jM 

P 


Ml ~ V 2 


sin (C — (5) + 2 cos 0i sin (0i + C + <5) 


(4.8) 


The above expression, except the finite term of the last line, presents certain divergences when evalu¬ 
ated at the end points of the integration interval, i.e. the functions log (|| 2 — 4>l) in the second line and 


arctanh 


l \/v 2 — v 2 


in the fifth line diverge at the upper and lower limits of the integration interval, respec¬ 
tively. In order to left with the finite value of the last line in the expression (4.8) we can choose the arbitrary 
parameters 6 , 0i in such a way that the coefficients of these divergent functions vanish identically. So, let 
us assume the following relationships between the parameters 


= 0, 


6 ucos (20i + £ + (5) + y^lV'ol |^4sin(£ + 0i) + 3 sin (0i + C + 2(5) - 3sin(0i - £)J 

Avy/p\ip 0 \cos (C + 0i) - 3vy/p\i /) 0 \cos (0i - C) + 3a/0#o| cos (C + 0i + 2(5) + 
3/31-001 2 sin (£ - 5) + 4/3|0o| 2 sin (£ + (5) — 6v 2 sin (20i + ( + 5) + 3/3|0o| 2 sin (20i + ( + 5) = 0. 

From the system (4.9)-(4.10) one gets the relationships 

0i = — ( — 5 — — arccos (2/3), 


( = arccos 


- A/i - (A 


V6 


(4.9) 

(4.10) 

(4.11) 

(4.12) 


Therefore, using the relations (4.9)-(4.10) in the expression (4.8) one gets a finite value for A f in (4.5) 
4|0o| 2 


u = — 


- \/vg — v 2 sin (C — (5) + 2 cos 0i sin (0i + a + <5) 


(4.13) 


where the parameters £ and 0i can be written in terms of v,v s ,5 through the relationships (4.11)-(4.12). 

Then, taking into account V = 2Mnls — §\//?| 0)| 3 from (3.26) and the above expression for A f written 
in terms of the parameters { 7 , 5 , v}, one gets the spatial shift from the equation (3.27) 


(4.14) 

(4.15) 


- 1 
LsXpert. — 

7 v 

'A( 7 ) cos 5 + B( 7 ) sin(5‘ 


A( 7 ) — — — 1 + 4 a/—4 + 57 2 
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( 4 . 16 ) 


B( 7 ) = - 2 v / 7 2 - 1 + V$\J- 4 + 57 2 + 2 v / 5 \/ 7 2 - 1 + 2\/5 

7 = = 7/#o|,> (4-17) 

where v s is the sound speed. Notice that the spatial shift Ax will depend on three independent parameters, 
{w,u s ,<5}. However, the parameter S can be fixed by considering the convergence of our series expansion to 
the value provided by the exact NLS phase shift in the limit v —> v s - 

On the other hand, the exact analytic expression for the spatial shift experienced by a black soliton after 
colliding with a gray soliton in the NLS model becomes (see (A.lO)-(A.ll)) 


a 1,7—1 

Ax analytic = X- log-— 

2 v s 7 +I 

= -L + oKrurt 

Whereas, the leading order term of the expression (4.14) in the limit v —t v s (7 >> 1) becomes 

(3/2) ( VI cos S + sin 5 


A.Xpert. — 


71 ? 




0\[lvY 


(4-18) 

(4.19) 


(4.20) 


Comparing the both expressions (4.19) and (4.20) one gets S = arccos y Then 

Ax nls = —— + 0[Vv)~ 2 }. (4.21) 

-yv 

Therefore, our expression (4.14) reproduces correctly the first order approximation in the parameter 
expansion 1 / (jv) of the spatial shift experienced by a black soliton after colliding with an shallow gray 
soliton in the NLS model. Let us evaluate the analytical results for some values of the parameters. The 
Fig. 1 shows a plot of the gray-black soliton collision in which the exact spatial shift using (4.18) becomes 
Axanaiy = —0.12437 which must be compared to the perturbation theory result Axnls ~ —0.1364, then 
this result provides an approximation within an error of 9%. This is acceptable since the velocity of the 
gray soliton being v = 4.64 (v s = 6, v < v s ) one has (vj) = 7.33, which is not expected to be a parameter 
characterizing an shallow soliton. The Fig. 4 shows the spatial shift Aa; vs (vj) for various values of 
the sound speed v s = 1, 6, 12. The continuous lines correspond to the exact analytical result (4.18) for 
v s = l(green),w s = 6(blue) and v s = 12 (red), respectively. The theoretical prediction (4.21) is plotted as 
a dashed line. The dashed and green lines exhibit an excellent agreement for (v'y) > 3. The inner plot 
shows that the expression in (4.21) exhibits a good convergence even for values (yv) ~ 3 provided that the 
sound speed satisfies v s < 1. Therefore, we can conclude that the leading order result (4.21) is in excellent 
agreement with the exact result of Aa; provided that v s < 1. Moreover, using the expression (4.17) of v s 
in terms of the parameters of the theory one can argue that a good convergence of the leading order result 
provided by (4.21) is restricted to a region in parameter space such that \tpo\ < l/VP- 

The behavior of the dark solitons in the sound speed limit (i.e. when the solitons become infinitely 
shallow, such that they reduce effectively to a pure plane wave) can be studied by considering it as being 
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equivalent to the Bogoliubov phonon, so the Bogoliubov-de Gennes equation (BdG) (3.12) associated in this 
case to the integrable NLS model can be considered as a scattering problem with the static black soliton 
potential. It is known that the dark soliton does not scatter phonons, i.e. the scattering of the phonon 
held with dark soliton potentials is rehectionless [7, 19]. This behavior has also been observed in BdG type 
equations with soliton potentials related to some integrable systems (see [20] and references therein). This 
feature implies that the spatial shift of the black soliton after colliding with an extremely shallow soliton 
must tend to zero in the sound speed limit. This behavior is consistent with our result above in this limit. 

The low velocity (v < NLS soliton collisions have recently been studied by considering the solitons 
as hard-splrere-like particles which interact through an effective repulsive potential [21]. The validity of the 
method has been checked by comparing to the exact analytical results for the collision-induced phase-shifts 
of the solitons, showing an excellent agreement in this regime. This analysis of soliton interactions has also 
been used to study a particular deformation of NLS model (2.1) ( /'[/] = \/l + 4 1 ) which considers an 
additional external harmonic potential [21]. The relevant theoretical results are in good agreement with 
the corresponding numerical simulations. In fact, at low velocities the solitons do not overlap completely 
during collision, instead they reflect from each other from the point of their closest proximity, so that they 
can be characterized by two individual density minima even at the collision region. However, high speed 
solitons (v > Hg-) are transmitted through each other, so they overlap completely in the interaction region 
and approach a single density minimum characterizing the location of the both solitons. This property has 
been used in our analysis above related to the regime in which the incoming shallow dark soliton possesses 
high speed and transmits through the stationary black soliton. 


4.2 Away from NLS: cubic-quintic non-linear Schrodinger (CQNLS) 

Consider the non-integrable cubic-quintic NLS model defined by 

i—^(a:,t) + - (>3\ip(x,t)\ 2 - ^0M)| 4 ) ip(x,t) = 0, (3 > 0,a > 0, (4.22) 

where we have considered /'[/] = (31 — j1 2 in (2.1). In the Appendix B we present the dark-soliton type 
solutions of this model, as well as the relevant expressions to be used in the equations (3.25)-(3.33) in order 
to compute the phase shift. So, let us compute the spatial shift experienced by a static black soliton (B.ll) 
after colliding with a shallow gray soliton (B.4). In the field space integration of Af (3.25) the field <j>i appears 
only in the expression of /'[|\l/o| 2 ] = /3|4/o| 2 — f 14/ 0 1 4 , where |4/o| 2 is given in (3.23); so using the relationship 
(B.l) and substituting the field 6 in terms of <f >2 from (B.8), the expression (3.25) can be written as 

M= _l fl [ m 02 #2 a23) 

v v a Jvu (IV'ol 2 - 4>2)V(4>2 - £i)(6 - 0f)’ 

where the integration in </>i in the interval [—IV’ol) iVbl] has already been performed and the integration limits 
of (j >2 is defined as the minimum and maximum intensities of the incoming dark soliton, i.e. [v^ft, |'0o|] (see 
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Appendix B). The functional /C becomes 


£[<h,C,0i,8\ = 


\HHHl - 2 /?) 

V& - Ci 


c 


l (a#o| 2 - 2/3) sin (£ + 6) - 


\H 




ja(38|?/>o| 4 


+ 


8O|?/>o| 2 02 + 15^) - 1O/?(4|'0 O | 2 + 3 ^ 2 ) + 10|V’o|[a|^o| 3 cos [2« + J]) + 
a\ipo\B 1 + 3a\ipo\B2 (3/3 — 4a|?/>o| 2 — 3a^2)2A 2 — 4aAi|^o| 2 cos (£ + 5) 


Ci - £2 

\/C 2 -Ci|V’o|sin(C + (5) | 


VC 2 - Ci 


C 


(4.24) 


M = \Ji 2 - \/Ci cos (01 - C) + ^I - Ci \/C2 sin (01 - C) 

A 2 = “ ^Iv^cos(0i -d) + ij<j>l -£i\/6sin(0i - S) 

Bi = IC2^ + CiWi - 2C2)] cos [2(0! - C)] - 2 y/(<% - C1XC2 - HV&fi sin [2(0! - Q] 

B 2 = + Ci( 0 l - 2C2)] cos [2(0! - 5)} - - C1XC2 - sin [2(0, - 5 )] 

c = -\M -Ci\/C2Cos (01 - C) + \/C2 - 0i-v/Ci sin (6»1 - £). 


Notice that in order to find the relationship (B. 8 ) between the fields </> 2 and 0 it is not necessary to have the 
explicit analytical solution of the gray soliton, since that relationship can be obtained from (2.12)-(2.14). It 
is clear that the integral in (4.23) diverges at the end points of the integration interval which are associated to 
the roots of the functional G[^> 2 ] (B.l). In order to extract a finite value for this integral, as in the computation 
of the usual NLS case above, one could examine the relevant divergent terms of the indefinite integral and 
choose convenient values for the set of parameters {£,5, 9{\ such that the definite integral becomes finite. 
However, in the present case this procedure turns out to be very cumbersome since it involves various types 
of Elliptic integrals and the relationship between the parameters would be complicated, therefore we will 
perform the computation considering a power series expansion in the perturbation parameter 1 /(^ 7 ). 

So, we make a relevant transformation in order to write (4.23) in the form (3.29)-(3.30). Taking the 
following transformation 


u = ($2 - £ 1 X 62 - <£!), 


(4.25) 


the expression (4.23) becomes 

AT 

{Hu]) 2 



p + + p 


1 + 


IC[H U . |,CA,<5] dy 


o [|V>o| 2 - {H u \) 2 ][p + - {Hu]) 2 ] 

3 v 1 \ 1 — u 2 11/2 


a{p+y 


T 


(4.26) 

(4.27) 


where p± = A |'0o1 2 A Remarkably, the overall factor l/(i> 7 ), which is the expansion parameter in 

the problem, appears in the expression (4.26). Then it will be sufficient to compute the zeroth order term of 
the integral in the power series expansion in the parameter 1 /(^ 7 ). However, this process must be performed 
carefully due to the divergence of the integrand at u — l(</> 2 = IV’ol) which appears explicitly in (4.26) and 
the indeterminate values which may appear when some terms of the functional /C[</> 2) C 0i, <5] vanishes at this 


18 



point. Notice that this divergent point is also achieved in the sound speed limit A- —> 0, £1 —> \i/jo\, £2 —> 
alJoP ’ ~^ l^°I ( see (B-9)). In fact, the functional 1C[(j) 2 , C, 9\,8\ in (4.24) contains some terms proportional 
1 ° \/02 — £1 which vanishes in the sound speed limit. 

Since we are interested in the sound speed limit of the integral in (4.26) let us write all the parameters 
as the series expansion x = Xq + ^ (x = {5, 9\, (, ^1 , C 2 }) and consider the series expansion of the 

integrand as H[<h[u],C, OiJ, l] = in powers of M 

^1)^)7] = 'H[<j> 2 [u = 1], Co, 9 10 ,8 0 ,7 ->• 00 ] + H'[<j) 2 [u = 1],Co,^io,<5o,7 ->• 00 ] (-) + ••• (4-28) 

7 

where H'lfalu = 1 ], Co, # 10 , <^o, 7 00] = 2W = 1], Co, ^ioA, 7 00], r = I/ 7 . 

We are interested in the zeroth order term of the series expansion above. For simplicity in the compu¬ 
tations we will choose some set of parameters {<5, 9 1 , £} in the sound speed limit (7 — > 00 ). In fact, a huge 
simplification is achieved for the choice 

c = — tt/2 + 0[7' 2 ], 6»1 = 0 + C>[ 7 - 2 ], 5 = 0 + O [ 1 ~ 2 }. ( 4 . 29 ) 


For the value 0\ « 0 the phase 9(z) of the gray soliton will remain as an odd function in the moving 
coordinate system, up to the lowest order (see the discussion below (B.6)). Moreover, the set of values (4.29) 
renders finite the lowest order term of the integration in (4.26). Considering the zeroth order contribution 
from (4.28) and performing the integration in (4.26), one has 


AT = - (■—) x (4.30) 

W ) 

_ 8q 3 /2 |V’o| 3 ^//3 - q|V>o| 2 (45 /3 2 - 110a/3|V> 0 | 2 + 63a 2 |^ 0 l 4 ) _ 

5 v / 3(3/3 - 4a|V’o| 2 ) 2 [v / «/3|V , o|v / 3v / / 3 ~ alV'ol 2 - (P ~ 2o#o| 2 )(3^ - 2a|^ 0 | 2 )tan~ 1 ( ^ 


From (B.13) one gets V = 2Mcq and therefore the spatial shift, using the relationship (3.27) becomes 


Axcqnls = 


J_\ _ 8A 3 / 2 y / T^A(45 - 110A + 63A 2 ) _ 

v "f) 5^(3 - 4A) 2 [(3 - 8A + 4A 2 ) arctan ( V3 ^_ A ) - v / 3\/ A ( 1 - A)] 
c(A) 


A = 


P 


v 7 


0 < A < 0.75. 


(4.31) 


This is the leading order in — expression for the spatial shift experienced by a black soliton after collision 


with the relevant gray soliton in the CQNLS model. Some comments are in order here. First, notice the 
appearance of the parameter A, which is defined as a relation between the quintic coupling a over the cubic 
coupling p constants multiplied by the squared amplitude |'0o| 2 - In this way up to the first order in 1 /(^ 7 ) 
the spatial shift Aa: does not depend on the individual parameters of the theory, rather on the effective 
parameter A which is a function of them. The Fig. 5 shows A xcqnls for various values of the effective 
parameter A. Second, the above spatial shift Aa; reverses its sign from negative to positive values at the 
point A = 0.654221 as it is shown in Fig. 8. This implies that there are some sets of values of the parameters 
where Aa: = 0. Third, in the limit A —> 0 one has Aa: = — (3 ^ 2 - ) (c(0) = 3/2), i.e. this is the spatial shift in the 
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sound speed limit when the parameters satisfy a|'i/ , o| 2 << /?• There is another value of A for which one has 
the same value of c(0), i.e. A = 0.526644 (the spatial shift for these values is plotted as a dashed line in Fig. 
5). Moreover, the expression c(A) above has a singular point at A = 0.75, then one must have (0 < A < 0.75). 
Therefore, our leading order result is meaningful in the regions of parameter space such that a < 4 |^| 2 ■ 
Fourth, the study of the relationship between the rate of convergence of the first order approximation in 
(4.31) and a set of parameters {cq/3, |V 7 o|} is interesting in its own right. The convergence check at a fixed 
(wy) can be performed by comparing (4.31) and its A dependence to the numerically calculated spatial shift. 
For a deformed NLS model this rate of convergence problem can be related to the study of the functional 
form of /[/] in eq. (2.1) and deserves a further exploration. 

Next, we present a brief discussion about the properties of the shallow CQNLS dark solitons. In ref. [14] 
it has been shown that the dark solitons become infinitely shallow at their bound for existence and that they 
are stable against weak perturbations. The problem of an infinitely shallow dark soliton and the outcome 
of its collision with the relevant black soliton can be discussed, as in the integrable NLS case above, in the 
context of a scattering problem of the Bogoliubov phonon and the static black soliton potential associated to 
the Bogoliubov-de Gennes equation (3.12) relevant to the CQNLS model. Recall that in the case of the BdG 
equation related to the integrable NLS model the phonon excitations are completely reflectionless against 
a dark soliton for any excitation energy. However, when a modified nonlinear term is present, the phonon 
acquires a nonzero reflection coefficient in general. This problem has been considered in [22] for the NLS 
model with general nonlinearity (2.1) and it has been shown that the property of perfect transmission of 
zero-energy phonon holds for a wide range of values of the velocity of the dark soliton which is used as 
the BdG equation potential. However, when a soliton velocity approaches a critical velocity this property 
disappears. In particular, a static black soliton does not scatter the zero-energy phonon. This feature implies 
that the spatial shift of the CQNLS black soliton after colliding with an infinitely shallow soliton must tend 
to zero. This behavior is consistent with our result above for the CQNLS model. 

5 Simulations 

The deformed NLS model (2.1) corresponds in general to a non-integrable model, so in the cases in which a 
dark soliton type solution exist in analytical form we do not expect to find an analytical expression for the 
collision of two solitons. That is the case for example in the cubic-quintic and certain cases of the NLS with 
saturable nonlinearities. However, the analytical expression of one solitary wave can still be used to provide 
an initial condition for the interaction of two solitary waves. So, in those cases we can take as the initial 
condition two one dark solitons some distance apart. 

After the fast-nroving gray soliton has moved away from the region occupied by the stationary black 
soliton, one is left with the stationary black soliton profile plus perturbations generated by the collision 
process. The spectrum of perturbations includes the spatial shift Aa; of the position of the black soliton 
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which we will compute numerically. The spatial shift can then be obtained by projecting the numerically 
calculated h(x,t) onto d x ipi(x) (and their relevant complex conjugates), so from (3.16) one can write 


A x(t) = 


if h(x,t)d x ipi(x) + h*(x,t)d x ip\{x) dx 


2 / d x ipi(x)d x ip$(x)dx 

where the stars stand for complex conjugation and h(x,t) can be obtained from (3.1) as 


(5.1) 


h(x,t) = ^{x,t) - ipi(x) - il> 2 {x,t) + ipo(x,t) (5.2) 

= ^(x,t)-Mx)e iwt e iC + \i> 0 \e iwt (e iS -e ie ~). (5.3) 

The function ^f(x,t) is provided by the numerical evolution of the full equation. We compute A x(t) in 
(5.1) after the fast moving gray soliton has already been collided and it is sufficiently far from the stationary 
black soliton, then the asymptotic value of the gray soliton has been taken as |'0o|e* 61 e lwt . It has been 
considered that the dark soliton is approaching from the negative direction of x, then after the collision one 
has the asymptotic value of the phase 0(—oo) = 0~. We consider a finite integration range in (5.1) in such 
a way that the derivatives d z \il)\\ and d z \il> 2 \ vanish outside the region of interaction of the solitons as in 
Figs. 1 and 2. Notice that the analytic result computed in (3.33) is time independent, whereas the numerical 
result (5.1) is not. There are a number of sources of time dependence. First, the finite box approximation 
for the numerical calculation renders the orthogonality of the zero-mode d x if>\ with the other vectors of the 
basis not to be exact. Second, the factors e lwt which appear in the dynamics of the individual dark solitons 
and the time dependence of the numerically calculated solution tE^a;, t). We extract the desired spatial shift 
by fitting the numerically calculated A x(t) to a straight line as 


A x(t) = Aa; + (A V)t, 


(5.4) 


where Ax will be the numerically estimated spatial shift of the black soliton, to be compared with the 
relevant analytical results in each case, i.e. expression (4.21) for the NLS case and (4.31) for the CQNLS 
case, respectively; AC represents the velocity change of the soliton. In general the last term in (5.4) accounts 
for the higher order [(fy) - ™, n > 2] time dependent terms which has not been considered in our analytical 
calculations. In the numerical simulations of the NLS and CQNLS models below we have verified that the 
velocity change behaves as (AC) = 0 + 0[(v^)~ 2 ]. In general the higher order terms will arise when the 
contributions from AW\h, A W^h and AWfh in the r.h.s. of (3.2) are properly taken into account, in addition 
to the source term S which has already been considered in our calculations. However, we expect that in the 
integrable NLS model AC = 0, to all orders in the perturbation expansion. The relevant developments for 
the higher order terms will require the knowledge of a complete set of orthonormal basis {^(a:)} in (3.11) 
which we postpone to a future publication. 

We performed extensive numerical simulations for the both the integrable NLS and the non-integrable 
cubic-quintic NLS models in order to check our results. The simulations are carried out by varying the 
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velocity of the incoming soliton while the other set of parameters of the model are maintained fixed, and 
afterwards by fixing the value (try) and varying the remaining constants of the model (e.g. {\ipo\,a, f3} in 
the CQNLS model). Next we present the numerical simulation results for the NLS model and the CQNLS 
models, respectively. 

NLS model. We have solved numerically the integrable NLS model defined in (4.1). We develop this 
procedure for an initial condition (t = 0 ) suitable for the interaction of gray and black solitons 


ip{x) 

C 


1 I {iv + P-2 tanh [P 2 (x + 72)]} , — L < x < x 0 

V 7 ? | Pi e tanh [Pi (a; + 71)], Xo < x < L 
v 


arctan 


LP 2 tanhP 2 (72 + xo)J ’ 


Pi cosh 


P2(72 + Xo) = p 2 cosh (Pi(71 + Xo)), 


(5.5) 

(5.6) 


where the solution (A.4) has been used, for the dark soliton one has v / 0 and for the black one v = 0. 
Despite the fact that one has an analytical solution for the gray-black system in the NLS model (A.5), which 
describes their collision at any time, the initial condition (5.5) considers the asymptotic form of each soliton 
when they are sufficiently far apart. The dark soliton is initially centered at —72 and moves to the right with 
velocity v, the stationary black soliton is centered at — 71 . Initially at the point Xo one has |^(xo)| = IV’ol 
and the both constant phases of the solitons match. The phase factor of the gray soliton e^ in (2.6) has 
been introduced in order to match the constant phase in (5.5). In the numerical simulation we will consider 
initially well-separated solitons, this amounts to considering the distance between them to be several times 
the width of the solitons, i.e. 1 71 — 72 1 > . Notice that the initial condition (5.5) must satisfy the 

boundary condition (2.3) at |x| —> L for |g/>o| = v ^ In the numerical simulation the domain 
considered is D = [—L,L] with L = 20, mesh size h = 0.022 and time step r = 0.0003. The domain D is 
chosen such that the effect of the extreme regions near the points x = ±L do not interfere the dynamics of 
the solitons. In our numerical simulations we have used the so-called time-splitting cosine pseudo-spectral 
finite difference (TSCP) method [23] suitable for the nvbc and the cw background at rest. The sound speed 
in the model ( 2 . 1 ) becomes v s = l^oI and the gray solitons of the model ( 2 . 1 ) must propagate with 
velocity |u| < v s . While we do not show the numerical results as in Fig. 2, we have also checked the results 
for various values of the sound velocity v s reproducing accurately the relevant figures. 
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Figure 4: (color online) spatial shift Ax as a function of (jv) in the integrable 
NLS model. The continuous lines correspond to the exact analytical result (4.18) 
for v s = 1 (green), v s = 6(blue) and v s = 12 (red), respectively. The theoretical 
prediction to first order of perturbation expansion (4.21) is plotted as a dashed 
line. Notice that for velocities close to the speed of sound (uy > 20) it shows a 
good agreement with the analytical result for the values v s > 1, whereas for low 
velocities it provides a good agreement for v s < 1 even for values vj ~ 3. 
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Cubic-quintic NLS (CQNLS). 

The numerical method TSCP is specialized in order to simulate the soliton collisions in the cubic-quintic 
NLS model (4.22). As the initial condition (t = 0) suitable for the interaction of gray and black solitons of 
the CQNLS model we consider the static solutions provided by (B.4)-(B.6) and (B.11)-(B.12), respectively, 
located far apart 


tp(x) 


A* 


gi+r 2 $2 tanh 2 [k 2 (x-x 2 )] iQ(x-x 2 ) P K 
l+r2 tanh 2 [k 2 (x—x 2 )] 

r+^tanh (x-xO] tailh fa (* ~ *1)] ^ e ' C ’ 


—L < x < xo 

Xq < X < L 


— arctan 



9{x 


X 2 ) = — arctan 



tanh [ k 2 {x — £ 2 )] j 


(5.7) 

(5.8) 


The gray soliton is initially centered at X 2 and moves to the right with velocity v = () 1 / 2 ! the stationary 
black soliton is initially centered at x\. Note that in the sound speed limit one has r 2 —> 0 (see eqs. (B.5) and 
(B.9)), so the constant phase above vanishes /r —> 0. Therefore it is necessary to introduce an overall factor 
e^ in the initial condition (5.7) in order to match the constant phase introduced for the black soliton (2.6). 
This parameter also appears in the perturbative computation of the spatial shift in (4.29). In the numerical 
simulation we will consider initially well-separated solitons, this amounts to considering the distance between 
them to be several times the width of the solitons, i.e. \x\ — x%\ > . Initially one has the conditions 

|^>(xo)| = \'4>o\i9 x \ij}(x)\ x — Xo = 0 and the both constant phases match at the point Xo■ Moreover, the initial 
condition satisfies the boundary condition (2.3) at \x\ —> L , i.e. |t/>(|ai| = L)\ = \ipo\- In the numerical 
simulation the domain considered is D = [—L,L\ with L = 30 (sometimes L = 20), mesh size h = 0.022 
and time step r = 0.0003. The domain D is chosen in order to give a well separated initial pulses and such 
that the effects of the extreme regions near the points x = ±L do not interfere the dynamics of the soliton 
interactions. 

The Figs. 6, 7 and 8 show the comparison of the analytical result (continuous lines) for the spatial shift 
as presented in (4.31) with the results obtained by numerical simulation (dots) of the gray-black soliton 
collision in the CQNLS model. The Fig. 6 shows the behavior of the phase shift as the velocity of the 
incoming gray soliton is varied. In this figure we used A = 0.2 which corresponds to the set of values 
\ipo\ = 2.75,/? = |V’o| 2 = 7.5625 and v s = 6.76411. Notice that the leading order result (4.31) reproduces 
the uq dependence of the spatial shift with a good approximation for high velocities given by (^ 7 ) > 20. 
Moreover, the difference between the numerically calculated spatial shift and the one calculated using (4.31) 
diminishes as (try) increases. The Fig. 7 shows how the analytic and numerical results approach each 
other as (vj) increases. In Fig. 8 we verify the A dependence of the spatial shift at a fixed value of 
(try) = 100. This value for ( 117 ) is fixed such that the analytical and numerical results have a good match for 
the particular set of values of the parameters chosen above. For the region X G h = [0,0.52] the numerical 
and analytic calculations are in good agreement. However, in order to achieve a good match in the region 
A € I 2 =< 0.52,0.75 > we choose another set of parameters, some values |t/»o| < 2.5, /3 = \4>o\ 2 have been 
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Figure 5: spatial shift Ax as a function of (71;) in the CQNLS model. The theoretical 
prediction based on the perturbation theory(4.31) is plotted for various values of the 
parameter A 0 = {0,0.526644}(c 0 = 1.5), A 2 = 0.409785(c 2 = 1.563), Ai = 0.62(ci = 
0.943). 



Figure 6: spatial shift Aa; as a function of (71;) in the CQNLS model. The theoretical 
prediction based on the perturbation theory result (4.31) is plotted as continuous 
line and the dots are the simulation results for the set of parameters a = 0.2, \tpo\ = 
2.75, P = |^o| 2 = 7.5625, v s = 6.76411, c(A = 0.2) = 1.53164. 
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Figure 7: We plot the product (v'j)Ax to show how the analytic and numerical 
calculations behave as (yu) increases in the CQNLS model. 



Figure 8 : We plot the product (uy) x A xcqnls as a function of A for (uy) = 100. 
The continuous line is the analytic result at leading order in l/(uy) and the dots 
correspond to numerical calculations. In the continuous line one has (—1.5) at A = 0 
and A = 0.526644, and the minimum value (—1.56319) is located at A = 0.409785. 
Notice that it changes sign at A = 0.654221. 
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used. Since the parameters we have chosen satisfy the relationship f3 = |" 00 1 2 one has A = - > A = a. 

So by varying A one is effectively probing the spatial shift behavior as if the quintic coupling constant a was 
varying. Therefore, the value of (vj) and the set of values of A at which the analytic result (4.31) provides 
a good agreement with the numerical calculations depend on the model parameters. Notice that the second 
interval I 2 contains the point A = 0.654221 G I 2 in which the analytic expression for the spatial shift changes 
sign (the interaction starts becoming attractive). Note also that at the limiting point A = 0.75 the leading 
order result diverges. 

We checked through our numerical simulations that the velocity change behaves as Av = 0 + 0[(v 7 )~ 2 ]- 
Finally, we conclude that the leading order calculation of the spatial shift is in good agreement with our 
numerical simulations for a wide range of parameter values. 

6 Discussions 

We developed a perturbative scheme to compute the spatial shift experienced by a black soliton after collision 
with an shallow dark soliton in deformed NLS models. The colliding shallow dark solitons pass through each 
other and for linearized perturbations the stationary soliton’s perturbations can be expanded as a power 
series in the parameter 1 /(^ 7 ). This is possible due to the fact that the overlapping region in space-time 
during collision possesses an area which is proportional to the factor A lnt ~ A-. This factor is extracted by 
simultaneous space-time and suitable field space transformations, which is provided for each deformed NLS 
model. We furnished a prescription to compute the leading order contribution to this free passage. 

We provided a closed expression for the leading order contribution to the spatial shift (3.33). This ex¬ 
pression depends on the the parameters of the theory, such as the coupling constants and the cwb amplitude, 
and the function f[I] which defines the equation of the deformed model (2.1). Remarkably, this result does 
not depend on the knowledge of the explicit solutions of the dark solitons. Two examples of applications are 
presented, the integrable NLS and the non-integrable CQNLS models. First, we have checked our analytical 
results with the integrable NLS model and found that they are in excellent agreement. Second, for the non- 
integrable CQNLS model we derived the leading order in ( 1 /( 117 )) contribution (4.31) which depends on an 
effective parameter A and checked with simulations. We found a good match with the numerical simulation 
results for a wide range of the parameter values. 

It must be emphasized the role played by the constant parameters associated to the phase factors of the 
black soliton and the cwb in the regularization procedure. This procedure has been performed in order to 
get finite values for the spatial shift in the both models presented above. 

Our analytical and numerical results show that qualitatively the shallow dark soliton collisions share 
some features with the kink-kink type collisions of the deformed SG model [11]. In fact, the absolute value 
of the spatial shift | Ax| decreases with increasing velocity in the both NLS and CQNLS models. In addition, 
in the CQNLS case and for fixed (vj) the spatial shift changes sign for certain value of the effective coupling 
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constant A. This resembles the behavior of the spatial shift in terms of the deformation parameter a for the 
deformed SG with potential V((j>) = (1 — cos^)(l — a sin 2 </>) [11, 12]. 

It would be interesting to see the effects of higher order terms and the remaining zero modes associated to 
phase, Galilean and scale invariance in future research. On the other hand, it is well known that the integrable 
models are characterized by the lack of velocity change after soliton collisions and that the quasi-integrable 
systems share some properties with their integrable counterparts. These properties of the quasi-integrable 
systems deserve future investigations in the framework of our formalism. 

For a non-integrable theory with solitary wave solution there appear some questions regarding the stability 
of the wave under small perturbations and the properties of the solitary wave collisions. Most of the known 
results concern some calculations in the regimes which are close to either integrable or related to a particular 
relationship between the parameters of the colliding solitons, e.g. high relative velocity, one of the solitons 
is significantly larger than the other one, fast thin solitons and slow broad solitons are among the cases 
considered in the literature. Related to these developments, when the wave configuration possesses a special 
space-time parity symmetry has been related to the concept of quasi-integrability [1]. In this context, the 
solitary wave interactions for the generalized KdV equation have been considered by Martel and Merle in a 
special regime [24]. We can also mention the work of G. Perelman on the bright solitary wave collisions of the 
deformed NLS model (2.1) with focusing nonlinearity, where one soliton is small with respect to the other 
[25] and the asymptotic approaches to describe the evolution and collision of three waves of the generalized 
KdV model by Omel’yanov [26]. 
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A Defocusing NLS: Exact analytic 2-dark soliton collision and 
spatial shift 

Here we present the one and two-dark solitons of the NLS model. The spatial shift experienced by a black 
soliton after collision with a gray soliton is computed. Consider the defocusing NLS equation 

id t ip + d xx ip - /3\ip\ 2 ip = 0, f}> 0. (A.l) 

The simplest solution of (A.l) has the form of a continuous wave 


ip = IV'ol exp [i(kx + wt + xq)], 


(A.2) 


provided that w = —At 2 /2 — (3\ipo\ 2 . In this work we will consider the continuous wave background at rest, 
i.e. k = 0. The system (A.l) is hyperbolic for small oscillations around \tp\ = \ipo\ and it has an associated 
sound speed given by v s = v^lV’ol [17]. This velocity plays an important role in the existence of solitary 
waves in the NLS model. 

The 1-dark soliton solution is given by [27] 


ip(x,t) 


1 

7$ 

i 

7 v 


t fl + y exp! zP(x-vt-xo)] 


}• 


\ l -|- expt 2P ( x v ^ 310 )] 
e lwt ^iv + Ptanh [P(x — vt — xq)]| 


y = e 2i \ 


?? = arctan {—P/v) 


(A.3) 

(A.4) 


This solution possesses three arbitrary real parameters, v,P and [3. The intensity function |^| moves at 
the velocity v, which is the velocity of the dark soliton. The dark soliton approaches constant amplitude 
\J ( P 2 + v 2 )/P as |x| —> ±oo. As x varies from —oo to ±00 the soliton acquires a phase 2 arctan (v/P) — n. 
Moreover, at the center of the soliton, i.e. for x — vt — x 0 = 0, one has that the intensity becomes \ip\center = 
vI\J~f3. This center intensity is lower than the asymptotic amplitude |'0o| and this property characterizes a 
dark soliton. Notice that this center intensity is controlled by the parameter 77 ; i.e. this parameter defines 
the “darkness” of the soliton. The soliton with v = 0 is the black soliton, whereas the moving soliton v ^ 0 
is the so called gray soliton. The speed of the soliton must satisfy |u| < v s . 

The 2-dark soliton solution is given by [27] 


ip(x,t) 
Ti (x,t) 


r 1 + 2/1 e ri ^’*) _|_ y 2 e I ' 2 ( x < t ') ± r yi 7/2 e ri ^ ,t ) e r 2 (a;,t) 


X -|- e riO,t) -|- _|_ J 'eri(x,t) e r 2 (x,t) 


2Pi(x — V\t — xi); T 2 (x,t) = 2P 2 (x 
(Pi - P 2) 2 + (Hi - V 2) 2 m . 

(Pi + P 2) 2 + (Vi — V 2) 2 ’ Vj 
— 7 t /2 < Sj < 7 t / 2 , j = 1 , 2 . 


-V 2 t 

, Sj = 


- x 2 ), 

arctan (—Pj/Vj), 


(A.5) 

(A.6) 

(A.7) 

(A.8) 


The parameters satisfy 


Pi + V 2 = P 2 + V 2 = (3\ipo\ 2 , w = -p\ipo\ 2 


(A.9) 
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The Vj are the soliton velocities and the P 3 are the parameters associated to the widths ( ~ ^p) of 
each soliton. Next, let us compute the spatial shift experienced by a black soliton after collision with a 
gray soliton. Let V\ = 0 be the black soliton velocity and V2 = v > 0 the gray soliton velocity. Assume 
Pi > 0 ,P 2 > 0 and that the gray soliton is initially located at x — > —00 then one has e T2 —>■ 0, therefore 
from (A.5) one can get ~ j }> which is the initial black soliton configuration at 

time ti- Consider the gray soliton is located at x —> +00 after it has already been collided with the black 

soliton, so e 1 " 2 —> 00 , and from (A. 5) one has ip(x,tf) ~ {—f —— ^ - P e 2 r 1 ^ |. Comparing this with 

the initial configuration one has that the spatial shift experienced by the black soliton becomes [28] 


Ax = 


1 7-I 

log— —r, 7 = 
2v a 7 +I 


'l-(v/v s )*’ 




(A.10) 

(A.ll) 


B Cubic-quintic NLS dark solitons 

In the following we will discuss some properties of the dark and black soliton type solutions of the cubic- 
quintic NLS model (CQNLS) defined in (4.22). Substituting (2.7) into the eq. (4.22) the corresponding 
expression in ( 2 . 12 ) becomes 

G[fa] = (B.l) 


Ci < l^ol 2 < £ 2 - 


Cl = 


C 2 = 


v 3 (j )2 V 

(3/3 - 2a\^o\ 2 ) - 1 /(3/3 - 2a|V>o | 2 ) 2 - 12 u 2 a 
2a ’ 

(3/3 - 2a|'0 o | 2 ) + sj(3j3 - 2q|V > 0 | 2 ) 2 - 12u 2 a 
2 a 


The set of parameters {Cij IV’ol) C 2 } are the roots of the functional G[<j) 2 ] and plays an important role in 
characterizing some properties of the dark soliton. Integration of (2.12) provides the dark soliton of the 
model [14] 


4%(z) = 


Ci + 7^2 tanh 2 [k 2 (z - z 0 )] 

1 + r 2 tanh 2 [k 2 (z — 2 0 )] 

IV’ol 2 - Ci , _ fa n — 


2 , fc 2 = V 3 V(C2 - |^0| 2 )(|^d| 2 - Cl)- 


Notice that k 2 characterizes the inverse soliton width and the root Ci is the minimum intensity (dip) of the 
dark soliton. The maximum intensity of the dark soliton approaches </> 2 (±oo) = |'0o|- Moreover, for v = 0 
the small root Ci vanishes and the dark soliton becomes a black soliton. 

Likewise using the solution (B.4) in (2.13) one has 


r 1 

9(z) = — arctan W — \Jr 2 tanh \ki(z — ^ 0 )] + #i, 

L V Ci 


where 9\ is a constant of integration. For 9\ = 0 the phase of the dark soliton approaches the constant values 
0 T = ± arctan as 2 —> + 00 , respectively, and it is an odd function 6(—(z — 20 )) = —9(z — zq) in the 
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moving frame of reference (z,t). Therefore, the dark soliton becomes parametrized as in (2.7) with (f> 2 {z) 
given by (B.4), 0(z) by (B.6) and 


W = -M o| 2 + ^IV’ol 4 


(B.7) 


In addition, the eq. (2.14) gives, upon integration and inversion the relationship 

1 + tan 2 (6 — 9\) 


= £2 


(B.8) 


' L|^ + tan 2 (0 — 0i) J 

The system (4.22) is hyperbolic for small oscillations around \ip\ = |'i/’o| and it has an associated sound 
speed given by v s = IV’ol v// 77 [IV’ol 2 ] = IV’oly//? — ajV’ol 2 [17]. Defining the parameter 7 = l/^/l — (^) 2 one 


can write the sound speed v —> v s (7 >> 1) limit of the parameters as 

7 ,. 2 

IV’ol 2 , 6^ 


2 • 


The stationary black soliton is obtained as the solution of the eq. (2.9). So, one has 

F [<t>i\ = yfd^ol 2 - 


4 ri tanh 2 (, k± z ) 

1 + rq tanh 2 ( k\ z) 

1 


_ IV’ol 2 


Z = —(3/3 — 2o;|g/>o| ), n = 

a 4 — 


ki = 


IV’ol 

Vs 


(B.9) 

(B.10) 

(B.ll) 

(B.12) 


The minimum intensity of the static black soliton <j> 1 is zero and it approaches |^o| as z —> ± 00 . Moreover, 
the eqs. (3.26)-(3.27) supplied with the eq. (B.10) give the energy of the black soliton of the CQNLS model 

1 


Mcq — 


4V 3 L 


IV’ol Vz - IV’o| 2 (C + 21V’o| 2 ) - Z(Z - 4|V’o| 2 ) arctan ( 




(B.13) 


Notice that this soliton mass formula reproduces the NLS black soliton mass (Mnls = |V/3|V’o| 3 ) in the 
limit a —> 0. 
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